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ABSTRACT 

In  a  report  ^  R.  M.  Lewis  has  shown  how  the  body  motion  and  geometric 
parameters  of  an  edge  located  on  a  torque- free  axi-symmetric  conducting  body 
could  be  determined  by  knowing  the  short  pulse  backscattered  cross  sections 
of  the  isolated  edge  at  two  different  polarizations  as  functions  of  time. 

The  purpose  of  this  report  is  to  indicate  another  more  direct  and  more 
practical  approach  to  the  inverse  edge  backscattering  problem  still  utilizing 
the  same  measured  data  as  Lewis's.  The  convenience  of  one  approach  over 
the  other  will  become  apparent  in  this  paper. 
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SECTION  I 


INTRODUCTION 


[1] 

In  a  report  R.  M.  Lewis  has  shown  how  the  body  motion  and 
geometric  parameters  of  an  edge  located  on  torque- free  axisymraetric 
conducting  bodies  could,  in  principle,  be  determined  by  knowing  the 
short-pulse  backscattered  cross  sections  of  the  isolated  edge  at  two 
different  polarizations  as  functions  of  time,  or  by  measuring  the 
time  history  of  the  short  pulse  polarization  scattering  matrix. 

The  required  short  pulse  must  be  short  enough  to  isolate  a 
particular  edge ,  but  not  so  short  that  in  effect  one  could  not  treat 
the  edge  scattering  as  essential  narrowband  scattering. 

Lewis  first  determines  the  local  direction  of  the  backscattered 
edge  by  maximizing  the  cross  section  as  a  function  of  the  polarization 
angle.  The  behavior  of  the  edge  direction  as  a  function  of  time  for 
torque-free  axisymmetric  conducting  bodies  can  then  be  utilized  to 
determine  the  kinematic  parameters  of  the  edge,  which  are  the  same  as 
those  of  the  body. 

Once  the  kinematic  parameters  of  the  body  are  known,  the 
geometric  parameters  of  the  edge  can  be  obtained  if  a  large  range  of 
aspect  angles  is  available. 

Lewis's  paper  mentions  that  there  are  many  different  ways  of 
determining  the  useful  parameters,  and  some  might  be  more  direct 
and  more  practical  than  the  one  introduced.  It  is  the  purpose  of 
this  paper  to  indicate  a  more  direct  and  practical  approach  for 
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determining  the  body  motion  parameters  of  torque-free  axisymmetric 
conducting  bodies  and  hence  the  geometric  parameters  of  the  observed 
edge  . 

The  method  utilizes  the  same  basic  data  as  Lewis's  method,  but 
the  theoretical  analysis  to  obtain  the  body  motion  parameters  is 
quite  different.  In  Section  II,  edge  backscattering  will  be  introduced 
in  the  light  of  our  approach  and  a  direct  way  of  attacking  the  inverse 
problem  will  be  considered. 

The  method  of  obtaining  body  motion  parameters  for  general 
torque-free  conducting  bodies  of  revolution  will  then  be  developed 
and  finally,  Section  IV  will  indicate  that ,  once  the  kinematic  parameters 
are  known,  the  method  for  determining  geometric  parameters  is  identical 
to  Lewis’s..  The  convenience  of  one  method  over  the  other  will  become 


apparent  in  this  paper. 


SECTION  II 


EDGE  BACK- SCATTERING 


r  2  31 

Using  geometrical  theory  of  diffraction  ,  one  can 
calculate^-  ^  the  diffraction  coefficients  for  electromagnetic 
backscattering  from  an  edge  to  be 
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qnd  where  the  angles  Y  and  0  and  the  subscripts  "1",  "2"  referring 
respectively  to  directions  ||  and  _L  to  the  local  edge  direction,  are 
shown  in  Figure  1. 

The  total  edge  backscattering  cross  section  for  incident 
polarization,  making  an  angle  0  with  the  edge  orientation,  is  given 
by 
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where  I  is  the  unit  vector  in  the  incident  field  direction.  N 

is  the  unit  normal  vector  to  the  edge,  and  p  is  the  local  radius 
of  curvature  of  the  edge  at  the  point  which  contributes  to  the 
scattering.  From  Equation  (2),  one  can  also  write: 
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The  b^  and  represent  the  measured  cross  sections  of  the  edge 

for  polarization  along  and  perpendicular  to  the  edge  direction; 

being  more  directly  measurable  than  the  diffraction  coefficients, 

they  will  be  used  throughout  this  paper. 

If  the  edge  were  a  circular  edge  on  a  body  of  revolution,  one 

sees  from  Equation  (4)  that  bx  and  b2  contain  information  about  the 
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dynamics  of  the  body  in  the  term 
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and  in  the  term 
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(because  of  0) ; 


however,  the  sum^  +  b^j  which  is 
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contains  information  about  the  dynamics  of  the  body  only  in  the  term 
1  ,  simplifying  considerably  the  inverse  problem. 
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Indeed,  one  can  write 
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b  +  b  =  g(Y  ,P)  /5TT 
VN-I 


(6) 


where 


(cos  TT/q  -  1) 


-1 


(7) 


is  a  function  of  the  geometrical  parameters  of  the  edge  and,  for  a 
circular  edge  it  is  a  constant. 

Using  a  short-pulse  radar,  it  is  possible  to  measure  the 
scattering  of  the  isolated  point  on  the  edge  contributing  to  the  back- 
scattering.  As  the  circular  edge  on  the  body  is  viewed  at  different 
directions,  the  point  which  contributes  to  the  backseat tering  moves 


A 


along  the  edge  and  hence  the  rotation  of  the  unit  vector  N  must 
properly  be  accounted  for. 

A 

A  mathematical  expression  for  N  can  easily  be  derived;  the 
unit  vector  must  lie  in  a  plane  formed  by  the  axis  of  symmetry  and 
the  line  of  sight  to  the  particular  edge,  and  thus  can  be  written  as: 


A 


(8) 


N  -  <^1  +  C2L 


where  L  is  the  unit  vector  along  the  axis  of  symmetry  of  the  body  of 


revolution  (Figure  2)  .  The  constant  C^,  and  C9  can  be  determined 
from  the  following  two  conditions: 
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Figure  2.  Geometrical  Relation  Between  the  Unit  Vectors  I 


The  first  condition  gives 
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and,  utilizing  the  second  condition,  one  obtains 
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therefore 
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The  inspection  of  Figure  2  sets  the  sign  at  plus  so  that 


N 


Using  the  above  expression  for  N  in  Equation  (6) ,  one  gets 
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and 


=  g2 (Y ,p)  ■■■■  . — 

:  *  *  9 
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(16) 


2 

Thus,  knowing  (b^  +  b^)  as  functions  of  time,  one  would  expect  that, 

if  the  body  of  revolution  were  undergoing  simple  torque-free  motion, 

the  above  expression  should  give  the  body-motion  parameters.  Before 

considering  a  method  of  obtaining  these  parameters,  let  us  investigate 

2 

first  how  from  the  measured  data  one  can  determine  (b^  +  b^)  as  a 

2  2 

function  of  time,  where  b^  and  a re  the  cross  sections  for  linear 

polarizations  along  and  perpendicular  to  the  edge  orientation. 

The  measurement  procedure,  which  we  establish,  is  the  following: 
not  knowing  the  edge  orientation,  the  radar  performs  a  two  step 
measurement  sequence  at  each  instant  of  time.  First,  it  transmits  a 
signal  which  is  linearly  polarized  along  a  horizontal  basis  vector 
and  measures  the  scattering  amplitudes  and  in  the  dual  linear 

polarization  receivers.  Second,  it  transmits  a  signal  which  is 
linearly  polarized  along  a  vertical  basis  vector  and  again  measures 
the  scattering  amplitudes  $22  linear  polarization 

receivers . 

Thus  at  each  instant  of  time,  the  radar  measures  the  polarization 
scattering  matrix  [s]  for  the  isolated  edge, 


Ls] 


11  S21 
12  S22 


(17) 
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Assuming  the  propagation  and  the  target  media  to  be  isotropic,  the 
scattering  matrix  will  then  be  symmetrical,  namely 


S12  S21 


(18) 


If  the  radar  had  the  horizontal  polarization  basis  vector 
aligned  along  the  edge  orientation,  then  according  to  the  geometrical 
theory  of  diffraction^,  the  scattering  matrix  [S']  would  be 
diagonal 


where  the  b's  are  the  quantities  which  we  seek  to  determine. 

Nevertheless,  having  measured  [S]  ,  one  can  always  find  [S'] 
and  the  angle  $  between  the  radar  horizontal  basis  vector  and  the 
edge  direction  by  solving  the  classical  eigenvalue  problem. 

In  our  case,  however,  since  one  does  not  need  to  know  0,  and 
b  and  b2  separately,  but  only  their  sum  (b.^  +  b2> ,  the  trace 
of  the  diagonal  matrix  [S']  ,  the  problem  is  considerably  simplified. 

The  measured  scattering  matrix  [S],  being  the  scattering  matrix 
for  the  edge  direction  oriented  at  an  angle  <t>  with  respect  to  the 
horizontal  polarization  basis  vector,  can  be  obtained  from  [S']  by 
a  [-<)>]  rotation  of  the  coordinates  of  [S']. 
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In  matrix  form^the  transformation  is 


[s]  =  [R]  ,  [S']  [R] 


-1 


(20) 


where  [R],  =  the  rotation  matrix 

<P 


COS  <j> 

sin 

-sin  <J> 

cos 

(21) 


Using  Equations  (19,  20,  and  21)  one  can  then  obtain  that 


ra  +  g  cos  2 4>  g  sin  2<J> 


[S] 


(22) 


g  sin  2i]i  o  -  g  cos  2< 


^1  +  b2  bl  b2 

where  a  =  - - -  and  g  =  - - - 


(23) 


Obtaining  the  trace  of  [S]  one  notes  that 


S11  +  S22  =  bl  +  b2 


(24) 


The  result  is  not  unexpected  since  it  is  well  known  that  the  trace  of 
a  matrix  remains  invariant  under  unitary  transformation;  such  as 
rotation,  for  instance. 

One  then  obtains  (b^  +  b^)  by  simply  measuring  the  diagonal 
elements  of  the  polarization  scattering  matrix  [S]  and  adding  them. 


11 


In  the  next  section  we  will  consider  how,  knowing  (b^  +  b^) 
as  function  of  time,  it  is  possible  to  determine  the  body  motion 
parameters  of  the  edge. 
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SECTION  III 


DETERMINATION  OF  BODY  MOTION  PARAMETERS 


As  previously  stated,  the  time  history  of  the  known  quantity 
(b^  4*  b^)  contains  information  about  the  dynamics  of  any  isolated 
edge  located  on  a  torque-free  axisymmetric  object.  From  Equation  (16), 
which  is 

(b  +  b  )2  -  g2(y,p)  1  <16> 

12  ^  [  *  A  2 

\jl  -  (I.L)Z 

where 

g(Y,P)  (17) 

q  k 


A  A 

one  sees  that  the  term  I-L  must  next  be  expressed  in  terms  of  the 
body  motion  parameters  of  the  torque-free  body. 

This  was  done  for  general  torque-free  axisymmetric  targets, 
by  the  author  in  another  report;  from  Figure  3,  for  the  case  that 

A  A 

R  »  (size  of  the  target)  ,  the  unit  vectors  I  and  L  are 


(30) 
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where  R  is  the  range  of  the  object  ,  Z  and  X  are  the  projections  of  R 
along  the  instantaneous  axis  of  precession  of  the  object  and  J. to  it, 
0*  is  the  semi-angle  of  precession,  U)  is  the  precessional  frequency, 
and  \|f  is  the  phase  angle  indicating  the  initial  position  of  the 
precessing  motion  relative  to  the  radar  frame. 

A  A 


Then  I.L  = 


-  [z  cos  0'  -I-  X  sin  0*  cos(uut  +  i|f)] 
R 


Thus,  one  clearly  sees  that  all  body  motion  parameters  are 

2 

included  in  the  quantity  (b^  4-  b^)  vs.  time,  and  the  extraction  of 
these  parameters  is  our  next  task. 


There  are  several  ways  of  proceeding: 


a)  One  way  of  obtaining  the  kinematic  parameters  is  to 


consider  the  quantity 
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which  in  terms  of  I-L  from  Equation  (15)  is: 
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Since 


(I-L)^  =  [ZcosB1  4-  Xsin0  '  cos  (Wt  4-  \|f)]‘ 
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The  known  quantity  S  is  then: 
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A  harmonic  analysis  will  give  the  precessional  frequency  and 

the  phase  angle  y.  Furthermore,  the  amplitudes  of  the  constant 

term  A  ,  of  the  1st  harmonic  An  and  the  2nd  harmonic  A  ,  which 
o  J-  A 

are  respectively 
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A,  = 


1  * /  s 
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R2 


05) 


A.  = 
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together  with  the  equation 


2  2  2 
R  =  Z  +  X 


(37) 


could  be  used  to  solve  for  the  four  unknowns : 


g  (y,p)  ,  Z,  X,  9'  . 


term,  it  could  be  exceedingly 


The  above  approach,  in  contrast  to  Lewis's,  does  not  require  a 

large  range  of  measurements,  making  it  very  suitable  for  rapid 

determination  of  the  dynamics  of  unstabilized  torque-free  objects; 

however,  because  of  the  — — 

(bj+bp' 

sensitive  to  errors,  requiring  very  exact  measurements  of  (b^  +  b^), 
b)  Another  method  of  extracting  the  body  motion  parameters 
from  Eq.  16  and  which  might  not  be  as  sensitive  to  errors  as  the 
previous  approach  is  the  following: 

Since  in  most  cases 
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one  could  use  the  binomial  expans 
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for  the  case  that  I»L  is  small  or  the  case  when  the  body's  axis  of 
precession  is  almost  1  to  the  line  of  sight  and  for  a  sufficiently 
small  angle  of  precession,  one  can  assume  that  (I*L)4  and  higher 
order  terms  can  be  neglected,  so  that 

(bx  +  b2)2  a  g2(V,P) 
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4R 


(39) 


2 

By  performing  a  harmonic  analysis  of  (b^  +  b^)  versus  time 
the  precessional  frequency  tl)  and  the  phase  angle  of  precession  could 
be  obtained.  Again,  by  knowing  the  amplitude  for  the  constant  part  Aq , 
the  1st  harmonic  A^,  and  the  second  harmonic  A^ ,  one  could  solve  the  set 
of  equations: 
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to  determine  the  remaining  unknown  body  motion  parameters. 

This  latter  approach,  although  not  as  sensitive  to  errors  as 

the  first  approach,  requires  a  more  specific  range  of  aspect  angles, 

A  a 

being  only  applicable  in  the  vicinity  of  the  region  where  I.L  =  0 

2 

or  where  (b^  4-  b^)  vs.  time  is  a  simple  periodic  function.  This 
again  limits  the  usefulness  of  the  approach,  when  trying  to  rapidly 
determine  the  dynamics  of  torque-free  bodies. 
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SECTION  IV 


CONCLUSION 


This  paper  has  illustrated  a  way  of  obtaining  body  motion 
parameters  different  from  that  indicated  by  Lewis^^;  however,  once 
the  body  motion  parameters  are  obtained,  the  method  of  determining 
the  geometrical  parameters  of  the  edge  are  identical  to  Lewis1 
In  short,  knowing  the  dynamics  of  the  edge,  one  utilizes  fully  the 
time  history  of  each  measured  cross  section  b^  and  b^  to  extract 
the  geometrical  parameters:  p ,  > ,  0  .  Indeed,  having  3  unknowns, 
it  is  sufficient  to  measure  one  of  the  b's  at  3  different  times  to 
obtain  a  set  of  equations  whose  solution  should  give  the  3  unknowns. 

It  is  very  clear  that  both  Lewis's  method  and  the  approach  described 
in  this  paper  are  very  useful  for  determining  the  kinematic  and 
geometric  parameters  of  torque- free  axisymmetric  targets;  furthermore, 
having  different  approaches,  the  radar  operator  has  more  leeway  to 
choose  the  one  which,  for  the  particular  case,  is  more  practical  and 
less  sensitive  to  errors. 
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